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Assertions that a new material may offer particularly advantageous properties should always be subjected to 
careful critical evaluation, especially when those properties can be affected by the presence of inclusions at 
trace level. This is particularly important for claims relating to new multiferroic compounds, which can 
easily be confounded by unobserved second phase magnetic inclusions. We demonstrate an original 
methodology for the detection, localization and quantification of second phase inclusions in thin Aurivillius 
type films. Additionally, we develop a dedicated statistical model and demonstrate its application to the 
analysis of BigT^ gFe! 5 2 Mno. 68 0 18 (B6TFMO) thin films, that makes it possible to put a high, defined 
confidence level (e.g. 99.5%) to the statement of £ new single phase multiferroic materials'. While our 
methodology has been specifically developed for magnetic inclusions, it can easily be adapted to any other 
material system that can be affected by low level inclusions. 



"A : 



bsence of Evidence is not Evidence of Absence" is a truism which confounds many analytical 
assertions. The confident statement of certainty: "There was nothing there" leads to the obvious 
. retort: "How do you know?". In many analyses, especially those using microscopes (optical or 
electron), it is only possible to look at restricted volumes, and the smaller the inclusion, or defect, or the lower 
the level of the impurity being sought, the smaller the volume that can reasonably be examined. The absence of 
any evidence for second phases, defects or impurities in these volumes can lead to the erroneous conclusion 
(perhaps driven by wishful thinking) that there are none in the sample as a whole, when in fact all that can be 
asserted with certainty was that there was nothing in the volumes examined at the limit of resolution for the 
technique used. And yet all is not lost, because if the detection thresholds of a particular method are well 
understood, the volumes examined are known and one makes some justifiable assumptions about the size 
distributions of second phase particles being sought, then one can apply statistical methods to determine the 
confidence level that surrounds such a statement. In other words, one can convert the assertion of "certainty-of- 
absence" to one which gives a level of probability that the density of the sought-for impurity or second phase 
inclusion is below a certain critical level. Conversely, it is possible to apply such methods to the design of 
experiments which will yield any required level of confidence. For example, it is possible to say that if one requires 
a confidence of 99.5% that the density of a certain impurity or defect of known composition and/or size is below a 
certain critical level in a specimen, then given the sensitivities of the analytical methodologies being used, the 
statistical analysis can be used to determine how large a total volume needs to be examined. 

If one considers different examples of materials analysis, there are different cases according to whether one is 
searching for impurity elements in a simple, single element material like silicon, trying to detect inclusions which 
are widely different in composition/structure from the parent matrix or looking for second phases that are similar 
in composition and/or structure to the parent. If the material is a pure element like silicon, for instance, claims 
with respect to a given level of purity can be verified relatively easily by using a host of different techniques such as 
SIMS (secondary ion mass spectrometry) and the level of impurities can be determined down to the level of ppb 1 " 4 . 
On the other hand, there are scanning methods which can locate inclusions and structural defects 5,6 , but, as 
commonly used, these are generally limited to detecting a minimum amount of 1-2 weight% of a particular 
inclusion type 7 . Atom probe tomography (APT) 8 and scanning transmission electron microscopy (STEM) 
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coupled with energy dispersive X-ray spectroscopy (EDS/EDX) 9 and/ 
or electron energy loss spectroscopy (EELS) 1011 provide techniques 
which can resolve single atoms spatially and elementally, but are 
limited to very small sample volumes (areas of a few nm 2 within a 
single atomic layer such as graphene film for example). Hence they 
cannot deal with a statistical distribution of small numbers of such 
entities dispersed over tens of urn 3 . There is a significant probability 
that the examination of a small volume of material at the high detail 
needed to see such very small inclusions may easily miss them, sim- 
ply because they are not present in the particular volume of material 
examined. In such a case it is prejudicial to say that such defects or 
inclusions are not present in the sample - one can only put a prob- 
ability upon them for their density being below a certain critical level. 

On the other hand, an aggressive search for new materials is in 
progress and the development of single-phase multiferroic (simulta- 
neously ferromagnetic and ferroelectric) oxides is a good example of 
this. These materials are of considerable interest for future memory 
and sensing applications 12 " 18 . Here, the assertion that a new material 
is a single-phase, room-temperature multiferroic may be confoun- 
ded by the presence of second-phase ferromagnetic inclusions that 
may have formed during the sample preparation. This is the single 
most critical factor in the analysis of such materials because the 
magnetic properties of the sample overall are highly susceptible to 
the presence of such inclusions. The difficulty in excluding second 
phase inclusions as being responsible for any observed ferromagnet- 
ism of a given sample is that some of the elements used to form the 
main phase could easily be responsible for formation of a ferromag- 
netic second phase. The differences between main and second phases 
could be in stoichiometry and/or crystal structure and both factors 
may have influence on the magnetic behaviour. 

This report demonstrates an original methodology for the detec- 
tion, localization and quantification of second phase inclusions in 
thin Aurivillius type films. Even when the compositional differences 
between inclusion and main phase are rather subtle, our methodo- 
logy provides a viable approach to not only identifying the inclusions 
but also to their localization and further in-depth structural 
investigation. 

Additionally, we develop a dedicated statistical model applied to 
the design of the analytical (EDX) measurements and demonstrate 
its application to the analysis of Bi 6 Ti 2 .8Fe 1 .5 2 Mn 0 68 O 18 (B6TFMO) 
thin films. We present a detailed account of the sample preparation, 
transmission electron microscopy (TEM) and energy dispersive 
X-ray analysis (EDX) that was necessary to put a high, defined con- 
fidence level (>99.5%) to the statement of 'new single phase multi- 
ferroic material'. 

While this methodology has been developed specifically to identify 
magnetic second phase inclusions, it can easily be applied to other 
material systems, given that the analytical measurement and statist- 
ical model is adapted to the special case, thus making the developed 
methodology comprehensively applicable. It closes an existing gap in 
the analysis of materials, between atomic scale analysis and the detec- 
tion of obvious large scale (micron size) inclusions. Due to the fast 
pace of development in the hard- and software for elemental detec- 
tion when using the electron beam as a probe, further development of 
the presented methodology can be sought by introducing necessary 
automation 19,20 . 

Statistical Analysis 

Our analysis requires a number of concepts from statistics, which we 
now introduce. We will start by analyzing the statistical significance 
of a single measurement of a given sample volume and at a fixed 
resolution, and then proceed to evaluate a series of measurements at 
varying volumes and resolutions. Finally we show how to conclude 
with a high, pre-set statistical confidence level (e.g. 99.5%) that a 
physical effect (which in the case at hand is the magnetic remanence 
of the sample) is not caused by undetected inclusions or defects. 



Single measurement. In general, it is impossible to prove with 
certainty that no inclusions exist in a given sample, unless the 
whole sample is scanned for inclusions. In order to nevertheless 
obtain statistically significant results, we use the method of 
statistical inference through the refutation of a null hypothesis. In 
the case at hand we assume that the level of concentration of a 
particular impurity exceeds a certain level, and then show that this 
assumption does not hold with a certain statistical level of 
confidence. More specifically, let us first assume that a particular 
type of second (or impurity) phase inclusion is present in the 
sample at a known (volume) density p and let us further assume 
that the individual inclusions do not interact with each other and 
are therefore independently distributed. We measure the number of 
inclusions found in a scan of a small volume V, as shown 
schematically in figure 1. The probability P(N, p) of finding 
precisely N inclusions is given by the Poisson distribution 

P(N,p)=^.e-" v . (1) 

In particular, the probability that no inclusions are present in a 
scanned volume V is given by 

P(0,p)=e- pv . 

If we now fix a desired confidence level y (in our case y = 99.5%) 
and define the value implicitly via 

£p(m,/£) = l-y (2) 

m = 0 

then it follows from (1) that for any value of the impurity density p 
bigger than p„, the probability of finding N or fewer inclusions within 
a volume V is below 1 — y = 0.5%. In the special case that no 
inclusions were found equation (2) can be solved explicitly and yields 

-/n(l-y) ^5.3 
Pu v V 9 

If a performed scan of the volume V yields N or fewer inclusions, 
we can therefore conclude that 

P<P N u (3) 

with a confidence of y = 99.5%. Subject to this confidence level, the 
p^ defined in equation (3) establishes an upper bound for the 
unknown inclusion density p. 

The grain size distribution function. So far we have assumed that 
there exists only one type of inclusion. In reality, inclusions can vary in 
grain size and composition. For simplicity, let us first study the case 
where the chemical composition of the inclusion is fixed (for example, 
Fe 3 0 4 as worst case scenario in our case) and all inclusions are sphe- 
rical. However, we do not make any prior assumptions with regards to 
the inclusion grain size. We can thus introduce a grain size 




Figure 1 | Schematic illustration of a scan of a volume V (small yellow 
box). What conclusions can we draw, if no inclusions (red balls) are 
detected during the scan? 
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distribution function p G {d), such that for small Ad the quantity 
p G (d)Ad is the density of inclusions with diameters in the range [d, 
d + Ad]. We can now generalize the upper bound (3) for the case of a 
size dependent distribution function as follows: assume that a scan of a 
certain volume V established that there are N inclusions present which 
are larger than a minimal detectable size d m . Then we can conclude 
that 

X 

dqp G (q)<Pu 



with confidence level y = 99.5%. 

Series of measurements. So far we have considered only one volume 
scan. In practice, any analysis will employ a number (K) of scans with 
different volume sizes V k possessing different minimal detectable 
grain sizes d k where the index k = 1, K enumerates the 
individual scans starting at large d k (i.e. d k > d k+ i) that yielded N k 
or fewer inclusions. Furthermore we assume that there is a physical 
upper bound d 0 > d x on the grain size. For a given grain size 
distribution function p G (d) we can define the density of inclusions 
bigger than d k via 



Pk= ^Pc(q)dq 



Then the probability of finding N k or fewer inclusions in one 
particular scan k is given by 

N k 

P k (N k )=J2 p (™,P k ) 

m = 0 

The probability of finding fewer than N k inclusions in a sequence 
of scans from k = 1, . . K is given by the product of the individual 
probabilities 

P = Pi(N 1 )P 2 (N 2 )---P K (N K ). 

Let us now assume that all individual scans establish that no inclu- 
sion was detected larger than the respective minimal grain size, d k . 
This is the case for the application at hand and therefore will now be 
discussed in more detail. The probability for this zero outcome is 
given by 

P(0)=P 1 (0)P 2 (0)---P K (0). 



The volume fraction (p k of inclusions between sizes d k and d k - Y can 
be calculated from the grain size distribution p G (d) assuming spher- 
ical inclusions 



<Pk = 



< 



dqp G (q) 



A- 



In this step we have employed the worst case scenario that the 
weight of the grain size distribution function is located at the large 
inclusion sizes. It therefore follows that there is an upper limit (p k>u on 
the volume fraction cp k of inclusions with diameter in the interval [d h 
d k -i] given by 



<Pk«Pk,u z 



-4- 



/«(l-y)«2.8 



4-i 



(4) 



which again holds with a confidence level of y = 99.5%. 

Here we have focused on spherical inclusions, but the generaliza- 
tion to other shapes is straight forward. For example, the assumption 
of cylindrical inclusions may be justified if the inclusion diameter 
exceeds the sample thickness. 



Effect of inclusions on remanent magnetisation. In the previous 
paragraph we established an upper limit of the volume fraction of 
inclusions in thin films. However, our aim is to establish an upper 
bound on the contribution of such inclusions on a measurable 
quantity. For multiferroic materials, the relevant quantity is the 
remanent magnetisation, which we will discuss now. Given the 
bound on the volume fraction, we can define the maximal 
contribution to the remanence from inclusions with diameters in 
the interval [d k , d k -{\ as 



M r Xu = <Pk,u M rA d k-l) 



(5) 



Here M r>F (d) is the remanent magnetization of pure Fe 3 0 4 grains of 
diameter d. For simplicity we assume that M r>F (d) vanishes for 
<5 nm. Furthermore we set M r)F (d) = 10 emu/g for 5 nm < d < 
20 nm and M r)F (d) = 20 emu/g for d > 20 nm which are 
conservative estimates. 

Finally we can establish an upper bound for the total contribution 
from inclusions to the remanent magnetisation as 



M nu -- 



max 

w u -,w K 



t ^ W k M rXu 

U=i 



with a confidence level y. Here, the parameters W l5 . . ., W K fulfill the 
constraint 



Therefore we can conclude that the contribution to the magnet- 
ization remanence from inclusions M r has an upper bound 

M r , M = max {MrXu} 

k=l...K 

with confidence y = 99.5%. 
Microstructural Analysis 

Guided by the statistical analysis detailed above, we designed our 
analysis measurements taking into account the sizes of the volumes 
scanned and the minimum size resolution of the EDX. Monte Carlo 
simulations 21 ' 22 of the electron beam matter interactions estimated a 
cone shaped interaction volume within the 200 nm film thickness 
with a base width of less than 200 nm and the EDX signal coming 
from a cylinder with roughly 15 nm diameter (interaction volume ~ 
35000 nm 3 ). The interaction volume is only dependent on the e- 
beam energy (20 kV), spot size (11.5 nm) and material (assuming 
that the main phase is the main contributor to calculated interac- 
tions) which all stay constant across all measurements. Taking all 
these factors into account a number of EDX maps were taken from 
the surface varying the surface area between 36 |im 2 and 
1,000,000 urn 2 (1 mm 2 ). As the mapping was done at a set resolution 
of 1024 X 886 pixels (max scanning resolution) we can determine the 
size of one pixel for each surface area. Hence the smallest detectable 
size of inclusion across all scanned areas is equal to 1 single pixel 
except for the smallest area where it is within 2X2 pixels. 

A set of 6 samples was made that went through microstructural 
analysis. All samples are Aurivillius type thin oxide films and their 
composition is summarised in Tab. 1. 

As described above, spinel-phase magnetite (Fe 3 0 4 ) 23 ' 24 is the 
worst case scenario for second phase inclusions in the type of samples 
we analyzed, but many other oxides exist that could be built from the 
elements of the main phase, and thus could occur as second phase 
inclusions. Some of these are also ferromagnetic. Hence it is neces- 
sary to scan for concentrations of all the magnetic elements that are 
present. The difficulty with these elements, however, is that they are 
also present in the main phase, and the ions in the main phase will 
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Table 1 


Summary of 6 samples 


that underwent microstructural 


analysis 






Sample 


Nominal composition 


Acronym 


1 


TbBi 5 Ti 3 Fe 2 Oi8 


TbBi5TFO 


2 


Tb 0 .5Bi5.5Ti3Fe 2 Oi 8 


Tb5Bi55TFO 


3 


Tbo.o5Bi5.95Ti3Fe20i 8 


Tb5B595TFO 


4 


Bi5Ti3Feo.7Coo.3O1 5 


BTF7C30 27 


5 


Bi 7 Ti3Fe 3 02i 


B7TFO 


6 


Bi 6 Ti 2 .8Fei.5 2 Mn 0 .68Oi8 


B6TFMO 



produce the majority of the X-ray signal detected for a given element. 
Figure 2 shows a typical result from a scanned area of 1600 urn 2 in 
B6TFMO (sample 6). At first, the elemental maps were colour coded 
(red - Ti, green - Fe or Mn, blue - Bi), added up into an RGB image 
and superimposed onto the SEM image, all done within the Oxford 
INCA software. Although an average recording of 5 X 10 9 counts 
over 72 hours improved the signal-to-noise-ratio by a factor of 71 
compared to a standard scan with an average of 1 X 10 6 total counts, 
and eliminated all the noise that was inhibiting the visibility of pos- 
sible inclusions in the maps achieved from shorter scans, the RGB 
colour coded images did not reveal inclusions instantly. The subtrac- 
tion of the Bi EDX signal from the EDX signals coming from the Fe 
and Mn allows us to detect subtle variations in the latter from the 
average value. The resulting maps are shown in figures 2d) Fe Ka 
map minus Bi La map and e) Mn Ka map minus Bi La map respect- 
ively which revealed small bright spots with higher Fe and Mn con- 
tent and larger (~1 um) pale areas with higher Fe content. The 
subtraction of the spectral imaging data was done under the assump- 
tion that the average Fe and Mn content in the main phase stays 
constant over the whole sample and all we are seeking to find are 
small changes from that average. Very important to note is the 
stability of the conditions during the 72 h scan. There is an approx. 
100 nm layer of carbon deposited over 72 h and this carbon is reg- 
ularly patterned. Each of the carbon dots represents exactly one pixel 
in the EDX measurement that was exposed approx. 1100 times. 

This method of subtracting the images to eliminate the signal from 
the Fe in the major phase was applied to all 6 samples. Except for the 
B6TFMO, all other analyzed samples showed potentially ferromag- 
netic inclusions of Fe x A y O z (see Figure 3) where A is substituted by 
Tb for the three samples containing Tb (samples 1-3), by Co for 
BTF7C30 (sample 4) and without A for B7TFO (sample 5). 

We also note that the inclusions are independently distributed. We 
have tested this by subdividing the images into partitions and count- 



ing the number of inclusions in each partition. The resulting count 
histograms, as represented by the typical example in Fig. 3f), show 
the expected Poisson distribution. It therefore follows that the spatial 
distribution of the inclusions is random and without detectable cor- 
relations, and justifies the basic assumption in the statistical analysis. 

The B6TFMO (sample 6) is interesting insofar as small bright 
spots (100-500 nm in diameter) with higher Fe and Mn content 
and larger paler areas (~1 urn in diameter) with higher Fe content 
have been found (Fig. 2, 4). Notably when comparing the Fe and Mn 
maps, the 'large pale areas' appear only in the Fe map (see Fig. 4). In 
order to investigate the microstructure and localization of some of 
these inclusions across the film thickness, and to correlate this to the 
magnetic response observed, a site- specific sample preparation was 
performed using a focused ion beam (FIB) (Fig. SI). At first the 
locations of a 'small bright spot' (labelled in yellow) and a 'large pale 
area' (labelled in blue) were marked with a pillar of e-beam deposited 
carbon (Fig. 5a). After the cross section lamella was successfully lifted 
out, attached to a copper half-grid and was polished to about 300 nm 
thickness the pillars were confirmed to be 250 nm in diameter and 
hence the location of the two areas was determined to this accuracy. 

With the two types of areas isolated, a STEM-EDX scan was per- 
formed as shown in figure 5b) and c). This measurement confirmed 
the absence of Mn excess in the 'large pale areas'. TEM imaging, in 
conjunction with selected area electron diffraction (SAED) and spot 
EDX measurements confirm this type of inclusion rather to be grains 
of the main Aurivillius phase with a 40% increased Fe content which 
has been associated previously to be the most probable cause of the 
magnetic behaviour of the sample 25 . 

TEM imaging, in conjunction with SAED and spot EDX 25 mea- 
surements confirm the 'small bright spots' as mangano-wiistite 
(Mno.56Feo.44O). This is reported to be non-ferromagnetic, antifer- 
romagnetic below 150 K 26 and hence not to be responsible for the 
magnetic behaviour of the sample. In conclusion, the inclusions 
found in the B6TFMO are not ferromagnetic. Hence their impact 
on the main phase magnetism is zero supported by statistical analysis 
which is elaborated in the Discussion. 

The link to the smallest volume in the statistical analysis for the 
B6TFMO sample consists of additional (168) bright field TEM 
images that were taken along a 30 um long lamella from a randomly 
chosen location for the structural analysis. There have been no inclu- 
sions observed in a volume of 0.32 um 3 . Every image shows smooth 
single crystalline Aurivillius phase material with grains touching 
each other smoothly and tightly (Fig S2). In our former studies of 
similar samples (BTF7C30 27 ) following such a "random" sampling 
approach, we were able to identify and image in detail the micro- 




subtract 



Bi La 





C 




Figure 2 | Example of EDX signal subtraction in case of B6TFMO sample: (a) Fe Ka intensity map, (b) Mn Ka intensity map, (c) Bi La intensity map. The 
results of subtraction of main phase bound Fe and Mn are shown in (d) and (e) respectively. 
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a) TbBi5TF0 b) Tb5Bi55TFO c) Tb5B595TFO 



d) BTF7C30 e) B7TF0 




!10 



Tb5Bi595TFO 



poisson 




0 12 3 4 
Inclusions per partition 



Figure 3 | (a)-(e) EDX surface maps for samples 1-5 respectively. Main phase bound Fe subtracted as per Fig. 2. The small bright spots visible are 
potentially ferromagnetic Fe x A y O z inclusions, (f) exemplary inclusion-count-rate histogram of (c) subdivided into 90 partitions, where the x-axis is the 
number of counts in a partition and the y-axis represents the number of partitions containing that number of counts. The theoretical poisson distribution 
is depicted for comparison. 



structure and composition of inclusions, similar to the type of second 
phase material being sought here. Hence for the sample under study 
in this report we can confirm that there are none. The resolution at 
the chosen magnification in the TEM was such that inclusions down 
to 2 nm would be identified, if present. 

Discussion 

The main question asked regarding the multiferroic oxide thin films 
is: 'Are there any ferromagnetic inclusions?'. This is the single most 
critical factor in the analysis of new single phase multiferroic materi- 
als because the magnetic response is particularly sensitive to the 
amount of such inclusions. To answer this question we developed 
a comprehensive framework based on linking the experimental EDX 
data from SEM and TEM scans with the numerical statistical analysis 
described above. The analytical framework we developed is schem- 
atically described in a flow chart (figure 6). The first step is in finding 
the optimum experimental analytical technique that can obtain spa- 
tially resolved maps of elemental compositions. EDX analysis is an 
obvious choice as it has the spatial resolution necessary to localize the 
impurities, and reasonably-good detection limits to most elements (Z 
> 10). Moreover, EDX analysis on SEM in the case of flat substrates 
does not require any extensive and elaborate sample preparation. 



The data acquisition (productivity of detection) has increasingly 
being improved by the introduction of large area solid phase detec- 
tors (SSD), increasing the solid angle of detection and automation 
procedures. In the next step, the obtained results of the experimental 
analysis can be treated as: YES/NO in regard to second phase 
detected. Once the choice is made the statistics need to be developed 
in order to supply the parameters that provide the input for the 
design of the experiment. In the section Statistical Analysis we dem- 
onstrate in detail the major steps in the development of the statistical 
analysis; starting at analyses of the statistical significance of a single 
measurement of a given sample volume and at a fixed resolution, we 
then proceeded to evaluate a series of measurements at varying 
volumes and resolutions. This methodology finally gives the possibil- 
ity to conclude, in this case with a statistical confidence level of 
99.5%, that a physical effect (e.g. the magnetic remanence of the 
sample) is not caused by undetected inclusions (NO choice feeds 
into the statistical analysis to calculate the upper limit of impact). 

In the case of inclusions being observed (YES outcome), the nature 
(structural, compositional properties) of these has to be determined 
in order to answer the next question in the flow chart: 'Do the 
observed inclusions have an effect on the magnetic remanence of 
the main phase?' If they do, the percentage of impact needs to be 



Fe b 



Q 3 



Small bright spot^ 

Q 



Q 4 

o° °« 

Q 3 G * Qi : 



O No large pale areas 

O Q 

_ Small bright spots'^ 

° ° G 



o 



Figure 4 | A comparison between Fe map (left) and Mn map (right) for the B6TFMO after subtraction of the Bi signal reveals the absence of the 'large 
pale areas' in the Mn map. 
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Small bright spot 




EDX map mix: 



Figure 5 | Elemental mapping of a cross section of B6TFMO prepared from the area shown in (a), (b) Fe is coded is green colour; (c) Mn is coded in green 
colour. 



calculated. For example, the presence of second phase magnetic 
inclusions at a volume fraction of 3.95% accounted for all of the 
observed magnetization in the BTF7C30 thin films reported prev- 
iously 27 . Though any trace of ferromagnetic inclusions within 
asserted multiferroic materials might raise concerns, the multiferroic 
property of a new material can still be established with a given con- 
fidence level, if the calculated impact of the inclusions seen is minor. 

In the case of no inclusions being observed (NO case) the contri- 
bution of inclusions to the remanence can be estimated statistically. 
As a worked example we have performed on the B6TFMO sample a 
number of measurements (labelled k = 1,...,8) at varying detection 
limits (<4 — 1 |im,... 2 nm) and scanning volumes (Vjt = 
200000 |im 3 ,...,0.32 (im 3 ) (see Tab. 2). It is not possible to scan the 
whole sample at the smallest resolution, and therefore the possibility 
for ferromagnetic inclusions cannot be excluded with certainty. 



Cost diagram 
Too coarse. 



^Optimum] 



Statistical 
upper limit 
of impact 




100%: Inclusions 
cause observed effect 



Figure 6 | The whole process flow of the new generic methodology for 
localization, identification and grain size distribution of inclusions in 
thin films with high confidence level. 



However, we can use the approach described in section Statistical 
Analysis to provide an upper bound for the contribution of undetec- 
ted inclusions to the remanence of the sample with a confidence level 
of 99.5%. The worst case scenario in the worked example of B6TFMO 
would be undetected Fe 3 0 4 inclusions and we therefore calculate the 
upper bound for this case. We assume that the remanent magnetiza- 
tion is directly proportional to the volume fraction of the inclusions. 
A direct correlation between the volume fraction of magnetic impur- 
ity inclusions and the observed magnetisation was previously 
demonstrated for the BTF7C30 sample 27 . We also assume that the 
magnetic inclusions do not interact with the Aurivillius phase matrix. 
This was justified by the lack of exchange-bias interactions at 2 K 
(±5 T) between the 4-layered BTF7C30 sample (Aurivillius phase is 



10 1 




io u 



0 50 100 150 200 250 
d [nm] 



10" J 



20 40 60 120 450 1000 3000 

d k [nm] 



Figure 7 | The blue bars indicate the maximal contribution M r ^ u to the 
remanent magnetization from possible unobserved inclusions with 
diameters in the interval [d^ d k - X ] as obtained from equation (5). The 
measured remanence of 22 memu/g for the B6TFMO sample is shown by 
the dashed red line. The inset shows the remanent magnetization M r>F of 
Fe 3 0 4 as function of grain size. Here the blue line shows literature data 
from reference 23 and reference 24. Because there are only a few literature 
values, we use a conveniently chosen (piecewise constant) green line which 
conservatively overestimates the remanent magnetization for all grain sizes 
d. 
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Table 2 Performed Vol 
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Scans for sample B6TFMO 
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antiferromagnetic at 2 K) and the ferrimagnetic CoFe 2 - x Ti x 0 4 
impurity inclusions (present at 3.95 vol.%) 28 . The effects of long- 
range magnetic interactions between unseen magnetic particles 
(<2 nm) were also excluded, since there was no evidence of a 
super-paramagnetic blocking temperature for B6TFMO 25 . 

To simplify the calculation we assume that the remanent magnet- 
ization of Fe 3 0 4 inclusions as a function of the grain size (M r>k>u ) is 
given by the green line in the inset of Fig. 7, which is chosen to be 
consistently above the literature values. We then calculate the values 
M n k,u which are the upper limits for the remanence contribution for 
each resolution interval \d h dk-i\ according to formula (5). The 
result is shown in the main panel of Fig. 7. We observe that all values 
M nk>u are well below the experimentally measured remanent magne- 
tisation of our sample (red dashed line of Fig. 7). For the upper limit 
of the remanence by unobserved Fe 3 0 4 inclusions we therefore find a 
value of M r>u = 0.69 memu/g with a confidence level 99.5%. This 
M r>u is only 3% of the measured magnetic remanence. This compares 
with a value of 13% as reported in 25 , for which the M r>u was calculated 
using a more restricted dataset as noted above. Following section 
Statistical Analysis this implies that the measured remanence is not 
due to undetected Fe 3 0 4 inclusions. This bound does not critically 
depend on the inclusions being spherical, for the errors made in the 
case of cubic, cylindric or ellipsoidal are minimal as the method of 
detection works with volume fractions independent of the shape of 
the inclusions. 

In conclusion we can say that we have developed a comprehensive 
and reliable method to detect, localize, isolate and analyse nm- and urn- 
inclusions in multiferroic thin films and calculated their impact on the 
properties of the main phase. More importantly this newly developed 
analytical methodology fills the gap between um inclusions and atomic 
scale inclusions. We are confident that this method can be applied to 
materials where a functional property (electrical, optical etc.), might be 
confounded by a second phase with dimensons in the nm and um 
range provided that the statistical design of the measurement is adapted 
to the particular problem and material system in question. 

Methods 

Thin film samples of Aurivillius type have been grown on c-plane sapphire substrates 
by chemical solution deposition method 25 ' 29 in order to investigate its potential as a 
single phase multiferroic material. 

Surface EDX mappings and cross section sample preparation required a thin Au 
coating (<20 nm) to reduce the high surface charging of the films. Cross sections 
were prepared using a DualBeam focused ion beam (FIB) FEI Helios NanoLab 600i 30 . 
After the Au sputter coating a 60 nm layer of amorphous carbon and 300 nm Pt layer 
were deposited within the DualBeam FIB by electron beam induced deposition and 
2 um thick Pt layer with ion beam induced deposition. These three layers have been 
grown for protection before the milling process. A 30 um long lamella was prepared 
and thinned down to 40-60 nm for the TEM analysis. The thinning at 30 kV was 
finished by polishing at 5 kV and 2 kV to reduce the ion-beam induced damage to a 
less than 2 nm thin layer on both sides 31 . 

TEM analysis was done in a JEOL JEM-2100 operating at 200 kV in bright field 
condition for imaging and selected area electron diffraction (SAED) mode to record 
diffraction pattern of single grains. Crystal Studio was used to determine the crystal 
structure by comparing the simulated electron diffraction pattern with the recorded. 



All EDX measurements were recorded in the FEI Helios NanoLab 600i operating at 
20 kV and 1.4 nA with an attached Oxford X-Max 80 detector. The spectra were 
analysed with the Inca software. 
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